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^ ■ Abstract 

H ; 

5t , We study tiie possibilities of spontaneous CP violation in the Next-to-Minimal 

Supersymmetric Standard Model with an extra singlet tadpole term in the scalar 
potential. We calculate the Higgs boson masses and couplings with radiative cor- 
rections including dominant two loop terms. We show that it is possible to satisfy 
the LEP constraints on the Higgs boson spectrum with non-trivial spontaneous CP 
violating phases. We also show that these phases could account for the observed 
value of ek- 



1 Introduction 



The understanding of CP violation, first observed in K decays [1], remains an open 
and most challenging question in particle physics. In the Standard Model (SM), CP 
violation arises from the presence of complex Yukawa couplings in the lagrangian. In 
electroweak interactions, CP violation originates in the misalignment of mass and charged- 
current interaction eigenstates, and is parametrized by the physical phase of the Cabibbo- 
Kobayashi-Maskawa (CKM) quark mixing matrix [2], 6ckm- Although SM predictions 
are in good agreement with experimental observations, including recent measurements 
at B-factorics, the SM amount of CP violation fails to account for the observed baryon 
asymmetry of the Universe [3]. Moreover, one is yet to find an answer for the strong 
CP problem, or in other words, to understand the smallness of the 9 parameter. Bounds 
from the electric dipole moment (EDM) of the neutron force this flavour-conserving CP 
violating phase to be as small as 10~^°, and this fine-tuning is most unnatural in the sense 
of 't Hooft [4] , since the lagrangian does not acquire any new symmetry in the limit where 
9 vanishes. 

In supersymmetric (SUSY) extensions of the SM there are additional sources of explicit 
CP violation, arising from complex soft SUSY breaking terms as well as from the complex 
SUSY conserving /i parameter. As pointed in [5], these phases can account for the values 
of the K and B meson CP violating observables, even in the absence of Sckm- However, 
the supersymmetric phases also generate large contributions to the EDMs of the electron, 
neutron and mercury atom. The non-observation of the EDMs imposes strong constraints 
on the SUSY phases, forcing them to be very small. Putting these new phases to zero 
is also not natural in the sense of 't Hooft. This is the so-caUed SUSY CP problem and 
many are the solutions that have been proposed to overcome it (for a review see Ref. [6]). 

An attractive approach to the SUSY CP problem is to impose CP invariance on 
the lagrangian, and spontaneously break it through complex vacuum expectation values 
(VEVs) for the Higgs scalar fields [7]. Thus, CP symmetry is restored at high energy and 
CP violating phases appear as dynamical variables. Spontaneous CP violation (SCPV) is 
also an appealing solution to the strong CP problem, since in this case one naturally has 
a vanishing 9 at tree level [8]. Further motivation to SCPV stems from string theories, 
where it has been shown that in string perturbation theory CP exists as a good symmetry 
that could be spontaneously broken [9] . 

SCPV requires at least two Higgs doublets. The Minimal Supersymmetric Standard 
Model (MSSM) is a very appealing example of a two Higgs doublet model. However, it 
is well known that, at tree level, SCPV does not occur in the MSSM [10]. On the other 
hand, radiative corrections can generate CP violating operators [11] but then, according 
to the Georgi-Pais theorem on radiatively broken global symmetries [12], one expects to 
have light Higgs states in the spectrum [13], which are excluded by LEP [14-16]. 

The case of the Next-to-Minimal Supersymmetric Standard Model (NMSSM) [17,18], 
where a singlet superfield is added to the Higgs sector, is more involved. In the usual 
Z3 invariant version of the model, where only dimensionless couplings are allowed in the 
superpotential, it has been shown that, although CP violating extrema are present at 
tree level, these are always maxima, not minima of the scalar potential, with negative 
squared masses for the Higgs states [19]. Hence SCPV is not feasible at tree level in the 
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NMSSM. Furthermore, as in the MSSM case, radiatively induced CP violating minima 
always bear light Higgs states which are difficult to accommodate with LEP data [20]. 
The possibility of SCPV has been studied in more general non-minimal models, where 
no Z3 symmetry is assumed and dimensionful, SUSY conserving terms are present in the 
superpotential [21-23]. In this case, it has been shown that SCPV is possible and could 
account for the observed value of Ek in the kaon system [21,23]. 

The first drawback of the general NMSSM with respect to the Z3 invariant version, 
is that it no longer provides a solution to the /i problem of the MSSM, which was one of 
the original motivations of the NMSSM. The second is that, in the absence of a global 
symmetry under which the singlet field is charged, divergent singlet tadpoles proportional 
to Mpianckj generated by non-renormalizable higher order interactions, can appear in the 
effective scalar potential [24,25]. Such tadpole terms would destabilize the hierarchy 
between the electroweak (EW) scale and the Planck scale. On the other hand, if Z3, or 
any other discrete symmetry, is imposed at the lagrangian level, it is spontaneously broken 
at the EW scale once the Higgs fields get non-vanishing VEVs, giving rise to disastrous 
cosmological domain walls [26]. It has recently been argued that using global discrete 
i?-symmetries for the complete theory - including non-renormalizable interactions - one 
could construct a Z3 invariant renormalizable superpotential and generate a Z3 breaking 
non-divergent singlet tadpole term in the scalar potential [27,28]. These models are free of 
both stability and domain wall problems, and all the dimensionful parameters, including 
the singlet tadpole, are generated through the soft SUSY breaking terms. 

The aim of this paper is to study the possibility of SCPV in the NMSSM with an extra 
singlet tadpole term in the effective potential, taking into account the latest experimental 
constraints on Higgs boson and sparticle masses, as well as the observed value of ex- In 
particular, contrary to what was asserted in Ref. [29], we obtain that in the Z3 symmetric 
limit, SCPV cannot be accommodated with the LEP exclusion limit on a light Higgs 
boson. On the other hand, we show that if Z3 is broken by a non-zero singlet tadpole, it 
is possible to spontaneously break CP, satisfy the LEP constraints on the Higgs boson mass 
and have Ek compatible with the experimental value. The paper is organized as follows: 
in section 2 we define the model and derive the Higgs boson mass matrix. The procedure 
to scan the parameter space and the resulting Higgs boson spectrum are discussed in 
section 3. Section 4 is devoted to the calculation ofEx- Finally, we present the conclusions 
in section 5. 

2 Overview of the model 
2.1 Higgs scalar potential 

In addition to the Yukawa couplings for the quarks and leptons (as in the MSSM), the 
superpotential of the NMSSM is defined by 

W^Higgs = XH1H2N + ^kN^ , (2.1) 

where Hi = {H^, H^ ) is the Higgs doublet superfield coupled to the down- type fermions, 
H2 — {H2, H2) coupled to the up-type ones, and TV is a singlet. Once EW sym- 
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metry is broken, the scalar component of N acquires a VEV, x — |(A^)|, thus generating 
an effective /i term 

H = Xx. (2.2) 

The superpotential in Eq. (2.1) is scale invariant, and the EW scale appears only 
through the soft SUSY breaking terms. It is also invariant under a global Z3 symmetry. 
The possible domain wall problem due to the spontaneous breaking of the Z3 symmetry 
at the EW scale is assumed to be solved by adding non-rcnormalizable interactions which 
break the Z3 symmetry without spoiling the quantum stability with unwanted divergent 
singlet tadpoles. This can be achieved by replacing the Z3 symmetry by a set of discrete R- 
symmetries, broken by the soft SUSY breaking terms [17,27]. At low energy, the additional 
non-renormalizable terms allowed by the i?- symmetries generate an extra linear term for 
the singlet in the effective potential, through tadpole loop diagrams 

T4adpoIe = -e'A^ + h.c. , (2.3) 

where ^ is of the order of the soft SUSY breaking terms (< 1 TeV). Since our approach 
is phenomenological, we take ^ as a free parameter, without considering the details of the 
non-renormalizable interactions that generate it. Likewise, we do not discard the singlet 
self-coupling term in Eq. (2.1) by imposing k — 0, which is possible once ^ 7^ [28], but 
rather assume k to be a free parameter. 

In addition to Vtadpoie, the tree level Higgs potential has the usual F and D terms as 
well as soft SUSY breaking terms: 

Vp = XmH,\'\N\^ + \H2\'\N\' + \H,\^\H2\') 
-A' {H^*H!^*H^H+ + h.c.) + K^\N\^ , 

Vd = ^-^^ - + ^-^\HlH,\\ (2.4) 

Koft = m]j^\Hi\'' + mjjjH^l^ + m%\N\^ 
+ {xAxNH^H^ + ^At^^iV^ + h.c. 

where gi and g2 are the U{1)y and SU{2)l coupling constants, respectively. In what 
follows, the soft SUSY breaking terms ■mHi,'mH2,'m'N, Ax, A,^ are taken as free parameters 
at the weak scale and no assumption is made on their value at the GUT scale. We assume 
that the lagrangian is CP invariant, which means that all the parameters appearing in 
Eqs. (2.3, 2.4) arc real. On the other hand, once the EW symmetry is spontaneously 
broken, the neutral Higgs fields acquire complex VEVs that spontaneously break CP. 
By gauge invariance, one can take — 0. The condition for a local minimum with 

{H2) = is equivalent to a positive square mass for the charged Higgs boson. The VEVs 
of the neutral Higgs fields have the general form 

{H^) = Vie'^P' , {Hi) = V2e''P'' , (A^) = xe'^^'' , (2.5) 
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where vi,V2,x are positive and ipi,ip2,ip3 are CP violating phases. However, only two of 
these phases are physical. They can be chosen as 



9 ^ (fii + (p2 + and S = 3(fi3 . (2.6) 
2.2 Minimization of the tree level potential 

Prom the tree level scalar potential in Eqs. (2.3, 2.4), one can derive the five minimization 
equations for the VEVs and phases vi,V2,x,9,5. They can be used to express the soft 
parameters mHniTiHi, f^Ni -^Xi in terms of Vi, V2, x, 9, S: 

^^ = m\ = -A2(x2 + v2sin2/5)-^M|cos2/3 

— Ax tan/? (kx cos(^^ — 5) + cos^) , 

^^ = ^ = -A2(x2 + ^;2cos2/3) + iM|cos2/? 

0V2 '2 

—\x cot /3 {kx cos{9 — S) + Ax cos 9) , 



dx 



m% = -AV-2«;V-A«;t;^sin2/?cos(^-(5) (2.7) 
■ sin2/3cos^ — kA^xcos^ H cos(5/3) , 



2x X 



dVtree „ . Kxsm(9 — 5) 

-^=0 ^ Ax = ^- 

o9 sm 9 

a\/tree _ q ^ _ BA/^xv^ sin 2(3 sm{9 -6) + 2^^ sm{5/3) 

dS 2kx'^ sin S 



with tan (3 — V2/V1, v — \/v\^ v\ = 174 GeV and Mz the Z boson mass. The above 
relations allow us to use tan /?, x, 9 and b instead of mij^, ■mH2,'>TT'N, Ax, A,^ as free param- 
eters. 

Once EW symmetry is spontaneously broken, we are left with five neutral Higgs states 
and a pair of charged Higgs states. The neutral Higgs fields can be rewritten in terms of 
CP eigenstates 



= e"^i + (5i + z sin pP) | , 
Hi = e"^' i^V2 + ^{S2 + i cos/?P)| , (2. 



N = e''^^|x + -^(X + iy)| , 

where Si,S2,X are the CP-even components, and P,Y arc the CP-odd components. 
Note that wc have rotated away the CP-odd would-be Goldstonc boson associated with 
the EW symmetry breaking. The mass matrix for the neutral Higgs bosons in the basis 
{Si, S2, X, P, Y) can be easily obtained: 
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, ,9 , ,9 9 ^ Xkx^ tan 3 sin 5 
Mfi = Ml cos^ /3 ^- 

Akx^ cot P sin 5 



= M|sin^/?- 



sin^ 



1 

MI3 = 4fi;2x2 + -Xkv^ sin 2/5sin(^ - 5) (cot ^ + 3 cot 5) + — (sin(5/3) cot 6 + cos(5/3)) 



2\k,x^ sin (5 



sin 2/5 sin 9 



Ml. = -A«^;2sin2/?f6cos(^-5)-9^-^V^(cos((^/3)-3sin((5/3)cot(5) 
,^9 /^^9 9 . Ara^sin^ 

(Sill (5 \ 
cos(^ - 5) + J (2.9) 



Ml^ = 

= 3A«;xi'sin/3sin(^ — 5) 

/ Sill (5 

Ml = 2A^x'u sin (5 + Arext; cos /3 cos (6' - S) A 

V sin 6* 

Ml^ = 

M25 = 3Afi;a;?; cos P sm{9 — S) 

= -AKa;'usin(^ - 5) 

MI5 = -2A«;'t;2sin2/3sin(^-(5) 

M^^ = Xkxv (cos{9 -S) + ) 
\ sm 9 J 

where we made use of the minimization conditions in Eq. (2.7) to ehminate the soft terms 
and simphfy the expressions. One can note that M^^ — M^^ = 0, which means that there 
is no CP violation in the Higgs doublet sector. On the other hand, CP violating mixings 
between the singlet and the doublets can appear, as long as ^ 7^ 5. It is easy to see from 
Eq. (2.9) that, if ,^ 7^ 0, one can evade the NMSSM no-go theorem for SCPV: it was 
shown in [19] that the tree level mass matrix for the neutral Higgs bosons always had one 
negative eigenvalue in the case of non-trivial phases of the VEVs. However, in our case 
the negative eigenvalue can be lifted up to a positive value provided ^ is large enough, 
due to the additional diagonal terms proportional to in MI3 and Mfg. Hence, SCPV is 
possible already at tree /ewe/ for 7^ 0. We will check this numerically in section 3, taking 
into account radiative corrections as well as experimental constraints on the Higgs boson 
masses. 



2.3 Radiative corrections 

It is well known that one loop radiative corrections can give large contributions to the 
Higgs boson masses in the MSSM [30] as well as in the NMSSM [31]. Furthermore, 
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they play a crucial role in the SCPV mechanism [11,20], as they generate CP violating 
operators. In what follows, we shall only consider radiative corrections due to top-stop 
loops. The one loop effective potential reads 

V.-,.., ^^,[Mt (,„g ^ - I) - 2Mt (log M _ I) I . (2,0, 

A4f and (i = 1, 2) are the field dependent top and stop squared masses respectively, 
and Q is the renormalization scale, at which all the parameters are evaluated. Q has to be 
of the order of the soft SUSY breaking parameters so that the tree level scalar potential 
has the supersymmetric form of Eq. (2.4). The field dependent stop mass matrix, in the 
basis (Ti,T^*), is given by 



ml + 1 i/o 1 2 ht {AtHi* + XH^N) 



ht {AtH^ + \H^*N* ) m^ + hl\Hl 



012 
2 I 



(2.11) 



with ht the top Yukawa coupling and rnqjUiT, At the soft terms for the stop sector. In 
the following, we assume mq — rriT = Msusy, as this choice maximizes the radiative 
corrections to the lightest Higgs boson mass. D terms are not taken into account since 
we do not consider radiative corrections proportional to the gauge couplings. At the 
minimum of the potential, the top-stop masses are given by 



2 7 2 2 • 2 ,o 

— h^v sm p , 

m| ^-^ = MsusY + T mtXt , (2.12) 

where 

Xt = ^/Al + A2x2 cot^ (3 + 2At\x cot (3 cos 9 (2.13) 

is the usual stop mixing parameter. This is not different from the CP conserving case, 
up to the phase 9 appearing in Xt- One can show that radiative corrections to the Higgs 
boson masses are minimized for = (minimal mixing scenario) and maximized for 
Xt = \/6MsusY (maximal mixing scenario). However, in our case it is not possible to 
take Xt — 0. In fact, the minimum of Xt is Xt — Xx cot /3| sin 6'| for At — —Xx cot /3 cos 9. 
On the other hand, Xt — -\/6Msusy is only possible if Ax cot P\ sm9\ < v^Msusy- If this 
is not the case, then the maximal mixing scenario is obtained for Xt — Ax cot f3\ sin9\ and 
At = —Ax cot f3 cos9 as in the minimum mixing case. 

The one loop terms give additional contributions to the minimization conditions of 
Eq. (2.7). In particular, since mj depends on the CP violating phase 9, the relation 
that gives Ax as a function of the VEVs and phases of the Higgs fields is no longer valid. 
Moreover, as noted before, the tree level minimization equations were used to simplify the 
Higgs boson mass matrix elements at tree level. Nevertheless, it is possible to keep the 
tree level mass matrix elements as in Eq. (2.9) and write all the one loop contributions 
as additional terms in the Higgs boson mass matrix 

Mi - + ^ SMf, . (2.14) 
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One then obtains 



5Mt, = -AVAf/K-,,m,-J 



2 2 
TTl- Tfl~ 

8Ml^ - _A2^2^(^ )+2AA2(7K^,m,-J + log^^ 

AVA^sin^^ ^, 

5M44 = r^- f{mi^,mi^) 

sm p 

^M^s = -AVcos2/3A2sin2^/K^,m,-J 
5Ml^ = -AxAAi A2 /(mj-j , J + AxAi ^((mij , m^J 

xx) cos 

^M^g = -A2x^;cos/3A2/(m,-^,m,-J + ^-^/iK^^^tJ (2-15) 

5M^5 = \^ XV COS 13 Atsin 9 f{mi^,mi^) 

5Ml^ = -A^; cos A2 /(mt^ , m^J + A^; cos /3Ai t/(mtj , m^-J 

.,^2 AxAfsin^s . Ax At sin ^ 

5M24 = /K-,m,-J-— -^,K-,m,J 



6MI^ — Xv cos 13 A^sm 6 /S.2 /(^ti^^ta) ~ Av cos /974( sin ^ g{mi^,mi^) 



SM^^ = AV cos' /5 A sin Ml /(mt^,mt J 

^,,9 A^xv cos sin^ 6* „, 

^Mi, = ^T^y /K-„m,-J 



where 



/(mi, 7712) 
g{mi,m2) 
/i(mi, 777,2) 



mo 



ml 



1 2 

7712 ~ ''^1 



log 



7772 
2 

777.7 



777; 



m 



2 log :i:2 ' 



7777 



2 2 
7772 ~ "^1 



1 ''^2 2 1 "^1 

log 7^ ~ "^1 log 7^ 



7717 



777.0 



Q 



1, 



(2.16) 



and 



Ai = Atcos^ + Axcot/3 , A2 = At cos 6* + Ax cot /9 cos 6* . (2.17) 

One can notice that, apart from a sub-leading term in SMf^^, all the scale dependence is 
hidden in the parameters scale dependence. The one loop terms of the neutral Higgs boson 
mass matrix in Eq.(2.15) differ from those given in Ref. [29], where one loop minimization 
conditions were not used to simplify the expressions. 
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Two loop corrections can also give substantial contributions to the Higgs boson masses 
[32,33] . Our analysis follows closely the results of [33] , to which we refer the reader for more 
details. Here, we consider the dominant two loop corrections which are proportional to 
ashf and h^, taking only the leading logarithms (LLs) into account. In this approximation, 
the two loop effective potential reads 

= 3 (^^) ' (32na. - I h^^ sin^ /? log' ^ . (2.18) 

One loop corrections to the tree level relations between bare parameters and physical 
observables, once reinserted in the one loop potential of Eq. (2.10), also appear as two 

loop effects. The dominant contributions are: 

(i) corrections to the kinetic terms of the Higgs bosons, which lead to a wave function 
renormalization factor Zjjo given by 

(n) Corrections to the top quark Yukawa coupling 

htimt) = htiQ) (l + ^ - 327r«,) log ^) . (2.20) 

The top quark running mass is then given by mt{mt) = ht{'mt)Z^JpV sin (3 and the rela- 
tion between the pole and running masses is, up to order m^°^^ = (1 + ^)mt{mt). 
Similarly, the relation between the stop pole and running masses is given by Tn?°^^ — 

(i + i^K-.K)- 

Once all these contributions are taken into account, one obtains a rather complicated 
5x5 mass matrix for the neutral Higgs fields, which can only be numerically diagonalized 
(cf. section 3). 

2.4 Charged Higgs boson mass 

Finally, we give the expression for the charged Higgs boson mass, which at tree level is 
given by 

mi± = - AV . 2.21 

^ ^ sin2/3sin^ ^ ^ 

The radiative corrections due to top-bottom and stop-sbottom loops read 

^ 9 3/i? XxAt cos 9 , , N r- 

Smjj^ = h{mi^,miJ + Sr , (2.22) 

where h{mi,m2) has been defined in Eq. (2.16) and Sr contains the scale independent 
terms 
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5r — 



1 

2 2~ 

?7i^ — mi 

t2 tl 



(m? 



log 



2 

OL 



+ 



t2 



2 2 

mi — m 

12 



167r^ m? + m? 

tl 12 



2 2 

mi 

tl t2 



log 



m 



mr 



(2.23) 



In the above, we have assumed that the bottom Yukawa coupling is small enough to be 
neglected, as justified by our choice of parameters, namely a low tan/3 regime. Therefore 
bottom squarks Bl and bu do not mix and m^^ = mg = Msusy- It is interesting to 
note that the term proportional to h{mf^,mfj in Eq. (2.22) exactly cancels the one loop 
correction to A\ due to the minimization equation of Vi_ioop as a function of 9. Hence, if 
one replaces Ax by a function of the VEVs and phases of the Higgs fields, one can simply 
rewrite the charged Higgs boson mass as 



w 



X'v' 



2Xkx^ sin 5 
sin 2/3 sin ^ 



tl t2 



(2.24) 



3 Mass spectrum 

In this section we investigate whether it is possible to have SCPV in the NMSSM with 
the extra tadpole term for the singlet, given the exclusion limits on the Higgs boson 
spectrum from LEP [14-16]. In order to do so, we perform a numerical scanning of the 
parameter space of the model. The parameters appearing in the tree level Higgs boson 
mass matrix are A, A;^, ^4^, m^/j, m//^! "^at, ■C- As seen in the previous section, we can 
use the minimization conditions of Eq. (2.7) to replace the soft SUSY breaking terms by 
the VEVs and phases of the Higgs fields. In the following, we will use the effective term 
defined in Eq. (2.2) as a free parameter instead of the singlet VEV x, so that the free 
parameters of the tree level Higgs boson mass matrix are now given by 

A , «; , tan/3 , , ^ , ^ , (5 . (3.1) 

Requiring the absence of Landau singularities for A and k below the GUT scale [Mqut ~ 
2 X 10^^ GeV) imposes upper bounds on these couplings at the weak scale, which depend on 
the value of /i*, i.e. of tan /3 and the top mass [18]. For tan /3 < 10 and = 175 GeV, 
one finds Amax ~ 0.65 and Kmax ~ 0.6. This also yields a lower bound for tan/3, namely 
tan/3>2.2. Regarding the radiative corrections, we take Msusy = 350 GeV and assume 
the maximal mixing scenario for the stops, i.e. Xt = max(v/6MsusYi A* cot /3| sin^^|). The 
relatively small value for Msusy will be justified in the following section when we address 
the computation of Ek- 

We have performed a numerical scanning on the free parameters, which were randomly 
chosen in the following intervals: 
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Figure 1: Mass of the lightest Higgs state as a function of ^, for a m; 



pole 



175 GeV, 



-^susY = 350 GeV and maximal stop mixing. The other parameters are randomly chosen 
as in Eq. (3.2). 



0.01 < A < 0.65 , 0.01 < K < 0.6 , 2.2 < tan /3 < 10 , 

100 GeV < < 500 GeV , < ^ < 1 TeV , (3.2) 

— TT < 9 < TT , — TT < 5 < TT . 

For each point, we computed the Higgs boson masses and couplings by diagonalizing 
the 5x5 Higgs boson mass matrix, which was calculated taking into account radiative 
corrections up to the dominant two loop terms, as described in the previous section. The 
five mass eigenstates are denoted by hi,i — 1..5 with masses rrii in increasing order. We 
also computed the charged Higgs boson mass, the stop and the chargino masses, and 
applied all the available experimental constraints on these particles from LEP [14-16, 34, 
35], as discussed below. 

The first important result is that we have succeeded in obtaining a large number of 
points that complied with all the imposed constraints for any values of the CP violating 
phases 9 and d, i.e. it is possible to have SCPV in the NMSSM. For this result to be valid, 
the presence of the extra tadpole term ^ is crucial. This is easily understood from the 
fact that in the hmit where ^ goes to zero, SCPV is no longer possible at tree level [19], 
and although viable when radiative corrections are included, the Georgi-Pais theorem [12] 
predicts the appearance of light Higgs states in the spectrum, already excluded by LEP. 
However, when ^ 7^ 0, CP can be spontaneously broken at tree level, as noted in the 
previous section. In Fig. 1 we display the mass of the lightest Higgs boson, mi, as a 
function of the tadpole parameter ^, with the other parameters randomly chosen as in 
Eq. (3.2), for a set of approximately 10^ points. We can see that small values of ^ are 
associated with a very hght mass for the hi state. Such light Higgs states are not excluded 
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Figure 2: Reduced coupling of the lightest Higgs state to the SM gauge bosons as a 
function of its mass for the same points as in Fig. 1. The solid line indicates the LEP 
exclusion limit [14]. 



by current experimental bounds as long as their reduced coupling to the gauge bosons 
is small enough. Indeed, the LEP exclusion limits on Higgs boson production gives an 
upper bound on the reduced coupling of a Higgs boson to the gauge bosons as a function 
of its mass. The reduced coupling Ri is defined as the coupling ZZhi, divided by the 
corresponding Standard Model coupling: 

Ri = hii cos [3 + hi2 sin (3 , (3.3) 

where hii,hi2 are the 81,82 components of the Higgs state hi, respectively. One has 
< l-Ril < 1 and unitarity implies 

E^' = l- (3.4) 

1=1 

In Fig. 2 we plot the square of the reduced coupling of the lightest Higgs state, Rf, 
versus mi for the same set of points as in Fig. 1. We also display the LEP exclusion 
curve [14], from which one can see that the light hi states are indeed not excluded. The 
presence of a light Higgs state with a small reduced coupling to the SM gauge bosons 
might prove difficult to detect at future colliders. However, it is worth stressing that 
in our results we always have at least one light visible Higgs state with a large reduced 
coupling to the gauge bosons (-R^>0.5) and a mass in the interval 112 GeV<mj< 150 GeV. 
The lower bound is fixed by the LEP limit as shown in Fig. 2. The upper bound can be 
understood from the following relation 

5 

Rjm] = cos^ /3Mi\ + 2 cos/5 sin (3M^2 + sin^ f3M^2 , (3-5) 

i=l 
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where Mf^, Mf^, M^^ are the Higgs boson mass matrix elements as given in the previ- 
ous section. The right hand side of Eq. (3.5) gives the usual NMSSM upper bound on 
the lightest Higgs boson mass [33], with the only difference being that the stop mixing 
parameter Xt now depends on the CP violating phase 9, as seen in Eq. (2.13). 

Assuming mf"^*^ = 175 GeV, Msusy — 350 GeV and maximal stop mixing, one obtains 
mi<125 GeV with this upper bound being saturated for small tan/?, namely tan (3 ~ 2.7 
[33]. One can check from Figs. 1 and 2 that this upper bound is only reached by a few 
points in our set. The dense band of points in Fig. 1 corresponds to cases where hi is a 
SM-like Higgs state with i?i ~ 1. In this case, the lower bound from LEP is mi > 114 GeV 
and the upper bound, from Eq. (3.5), is mi<125 GcV as explained above. 

We have also applied the exclusion limit from LEP on Higgs bosons associated pro- 
duction (e+e~ — > hA in the MSSM). This provides an upper bound on the Zhihj reduced 
coupling as a function of -|- mj for rrii ~ rrij [15]. The Zhihj reduced coupling is the 
equivalent of cos(/3 — a) in the MSSM and is here defined as 

R'ij — hi4{hj2 cos /9 — hji sin P) — {i ^ j) , (3.6) 

with /ij4 the CP-odd doublet P component of hi. We have also taken into account the LEP 
limit on the charged Higgs boson mass [16], which here simply reads m//± > 89.6 GcV. 
By inspection of Eq. (2.24), this limit translates into having opposite signs for the phases 
9 and 5. 

As we will see in the next section, charginos play an important role in the computation 
of Sk- The tree level chargino mass matrix in the (W , H) basis reads 

where M2 is the soft wino mass. We randomly scanned in the following interval for M2 

100 GeV < M2 < 250 GeV (3.8) 

and apphed the LEP bound on the chargino mass, Tn^±^ > 103.5 GeV [34]. Finally, 
we checked that for each point the stop masses also satisfied the associated LEP limit 
mf^ ^ > 100 GeV [35]. 

4 Indirect CP violation in mixing 

In the framework of the NMSSM with SCPV, all the SUSY parameters are real. Even so, 
the physical phases of the Higgs doublets and singlet appear in the scalar fermion, chargino 
and neutralino mass matrices, as well as in several interaction vertices. In this section we 
will explore whether or not these physical phases can account for the experimental value 
of Ek, Sk = (2.271 ± 0.017) x lO'^ pg]^ 

It is worth stressing that in this scenario the SM does not provide any contribution 
to the CP violating observables, since the CKM matrix is real. This can be clarified by 
noting that since N is a, singlet field, it does not couple to the quarks, and although both 
Higgs doublets do couple, the phase associated with these couplings can be rotated away 
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Figure 3: Box diagrams associated with the leading chargino contribution to Ek- 



by means of a redefinition of the right-handed quark fields. Since charged currents are 
purely left-handed, these phases do not show up in the Vckmj which is thus real. 

Let us now proceed to compute the contributions to the indirect CP violation param- 
eter of the kaon sector, which is defined as 

- • (4.1) 

ArriK ^ ' 

In the latter ArriK is the long- and short-lived kaon mass difference, and Aiu is the 
off-diagonal element of the neutral kaon mass matrix, which is related to the effective 
hamiltonian that governs AS — 2 transitions as 



M 
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2m 



with 



K 



o/A5=2 
^eff 



(4.2) 



Here q are the Wilson coefficients and Oi the local operators. In the presence of SUSY 
contributions, the Wilson coefficients can be decomposed as Ci = cf^ + cf +cf +cf + cf. 
As discussed in [23] , in the present class of models where there are no contributions from 
the SM, the chargino mediated box diagrams give the leading supersymmetric contribu- 
tion, and the AS* = 2 transition is largely dominated by the {V — A) four fermion operator 
Oi- Working in the weak basis for the W — H, rather than in the physical chargino basis, 
and using the mass insertion approximation for the internal squarks, we have verified that 
the Sk receives the leading contribution from the box diagrams depicted in Fig. 3. In the 
limit of degenerate masses for the left-handed up-squarks, Im AI12 is given by [23] 



Im A^i2 



{V*aVts)m^ \e'^ - cot/3 m^l 



37r2(mg)8 

X [^Au sm[^^ - e] (M?)^^ I{rw^rH,ri,^,rijj 



(4.3) 
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Figure 4: Values of sk as function of tan (3 (shaded area) for = 175 GeV, MgusY = 
350 GeV and maximal stop mixing. The other parameters are randomly chosen as in 
Eqs. (3.2, 3.8). 



In the above equation, fx is the kaon decay constant and mx the kaon mass [36]; Vij are 
the VcKM elements, whose numerical values (Vtd = 0.0066 and Vts = —0.04) reflect the 
fact that we are dealing with a fiat unitarity triangle; {nig) is the average squark mass, 
which we take equal to Msusy; ~ wino mass, = fi is the higgsino 

mass and ip^ = arg(e*^m^y — cot/? i^hY- The non-universality in the LL soft breaking 
masses is parametrized by (M?) , which we choose taking into account the bounds from 



the analysis in Ref. [37]. As discussed in Ref. [23], a sizable non-universality in the soft 
trilinear terms is crucial, and is here parametrized by /S.Ajj = — ^'u- ■'■^ following, 
we assume AAu = 500 GeV. Finally, I is the loop function, with = mf/{mq)'^ [23]. 

By scanning the parameter space of the model, we have verified that one can find 
sets of parameters that satisfy the minimization conditions of the Higgs potential, have 
an associated Higgs boson spectrum compatible with LEP searches and still succeed in 
generating the observed value of Sk- From Eq. (4.3), it appears that Sk depends on 
1/MgusY- It is therefore difficult to saturate the experimental value of Ek for large values 
of Msusy- On the other hand, too small values of MgusY naight generate a light Higgs 
boson spectrum, already excluded by LEP. In order to accommodate both constraints, we 
took Msusy = 350 GeV, as already referred to in section 3. All other parameters are as 
in Eqs. (3.2, 3.8) and we assumed = 175 GeV and maximal stop mixing, as in the 

previous section. 

In Fig. 4, we plot the possible values of Ek as a function of tan (3. The maximal values 
of Ek are obtained for the low tan (3 regime, and the experimental bounds on ek require 
tan/3<3.8. Recall that the lower bound for tan /? is consistent with the analysis conducted 
in section 3. Apart from the explicit dependence of Eq. (4.3) on tan/?, one should bear 
in mind that there are also implicit dependences associated with the other parameters 

^Our conventions differ from Ref. [23] in that /ie// — 
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Figure 5: Contours for the maximum value of Ek in the X-k plane for mf°^ = 175 GeV, 
^susY = 350 GeV, maximal stop mixing and the other parameters as in Eqs. (3.2, 3.8). 



as well as with the various experimental bounds imposed on the mass spectrum. It is 
therefore difficult to reproduce analytically the observed upper bound on ek as a function 
of tan (3. 

In Fig. 5, we plot the maximal value of Ek in distinct regions of the \ — k plane. 
The remaining parameters are chosen in order to maximize Ek and still comply with the 
experimental bounds. As one can see from this figure, having ~ 2 x 10^^ is associated 
with values of k and A in the range [0, 0.6]. In other words, one can easily saturate Ek in 
a vast region of the parameter space. 

As expected from the inspection of Eq. (4.3), there is a strong dependence of Ek on 
the phases associated with the Higgs fields. In Fig. 6, we display contour plots for the 
maximal values of Ek in the plane generated by the phases 6 and 5. Although, a priori, 
all the values for the phases 9 and 5 in [— tt, tt] are allowed, it is clear from Fig. 6 that the 
saturation of the experimental value of ek can only be achieved for significant values of 
the singlet and doublet phases. This feature becomes more evident in Fig. 7, where we 
show the values of Ek as a function of the singlet phase, ip^, = S/3. The saturation of the 
experimental value of Ek requires the singlet phase to be \(ps\ > 0.15. We will discuss in 
the conclusions the implications of such large CP violating phases. 

5 Conclusions and discussion 

We have shown that in the framework of the NMSSM with an extra tadpole term for 
the singlet in the scalar potential, it is possible to have spontaneous CP violation while 
complying with the LEP constraints on the Higgs boson mass. Although in this scenario 
VcKM is real, the experimental value of Ek can be saturated in large parts of the parameter 
space of the model, without requiring any fine tuning. This may be of special interest for 
electroweak baryogenesis, as it has been shown that a strong first order electroweak phase 
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Figure 6: Contours for the maximum value of Ek in the 6-5 plane. All parameters as in 
Fiff. 5. 
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Figure 7: Values of Ek as function of the singlet phase (ps = 6/3 (grey area). All param- 
eters as in Fig. 5. 
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transition is possible in the NMSSM, as opposed to the MSSM, due to the additional 
trilinear Higgs boson couplings in the lagrangian [39]. 

The presence of the extra tadpole term ^ is crucial to concihate spontaneous CP 
violation with the LEP constraints on the Higgs boson mass. Moreover, the ^ term 
also solves the domain wall problem associated with the spontaneous breaking of the Z3 
symmetry, which is present if ^ = 0. On the other hand, since CP is spontaneously broken 
at the electroweak scale, CP domain walls may appear, which are cosmologically excluded. 
An elegant solution to this problem, along the hues of the Z3 domain walls solution, is to 
assume that gravitational interactions also explicitly break CP. In the low energy scalar 
potential, one then obtains an explicit CP violating tadpole term, i.e. a complex term. 
A phase in the ^ term would not change the results derived here, since its only effect is 
to generate a shift in the singlet phase 6. 

Regarding the large phase regime favoured by the saturation of ex, one should recall 
that 9 and 5 are flavour-conserving phases, and might generate sizable contributions to 
the electron, neutron and mercury atom EDMs. Although we will not address the EDM 
problem here, a few remarks are in order: first, let us notice that in the presence of a 
small singlet coupling A, as allowed in our results (see Fig. 5), the EDM constraints on 6 
become less stringent [38]. In addition, there are several possible ways to evade the EDM 
problem, namely reinforcing the non- universality on the trilinear terms (i.e. requiring the 
diagonal terms to be much smaller than the off-diagonal ones or having matrix-factorizable 
A terms), the existence of cancellations between the several SUSY contributions, and the 
suppression of the EDMs by a heavy SUSY spectrum [6]. In view of the considerably large 
parameter space allowed in our results, none of these possibilities should be disregarded. 

Concerning the other CP- violating observables, namely e'/e and the CP asymmetry of 
the Bd meson decay {aj/^Ks)-, it has been pointed out that this class of models can generate 
sizable contributions, although saturating the experimental values generally favours a 
regime of large phases and maximal LR squark mixing [21,40]. A complete analysis of 
these issues, as well as of the EDM problem in our model has yet to be done. 

In this work, we pointed out that the NMSSM with an extra tadpole term in the 
scalar potential appears to be an excellent candidate for a scenario of SCPV. Having a Z3 
invariant superpotential preserves the original motivation of the NMSSM to solve the 
problem of supersymmetry. The extra tadpole term for the singlet cures the domain wall 
problem of the Z3 invariant model and allows the spontaneous breaking of CP. In this 
model, one can simultaneously saturate Sk and obtain a particle spectrum compatible 
with the current experimental bounds. 

The presently available experimental data on CP violation does not allow to distinguish 
whether the CP symmetry is spontaneously or explicitly broken. With the advent of the 
LHC and eventually of the Linear CoUider, direct searches of the Higgs bosons, and the 
measurement the associated couplings to the SM gauge bosons, will allow to disentangle 
this model from the MSSM. 
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